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ABSTRACT 

We reanalyze our recently proposed mass matrix model based on 
spontaneously broken discrete chiral family symmetry taking into account 
the additional flavor changing neutral current constraint implied by the 
bound on the D 1 — D 2 mass difference, and including several corrections 
to our earlier analysis. When combined, the K\ — K 2 and D\ — D 2 con- 
straints force the masses of the Higgs particles that contribute most strongly 
to flavor changing neutral currents (the Higgs states) to lie above 17 TeV, 
well beyond the limit of validity of conventional perturbative Higgs physics. 
The analogous constraints on the masses of the rj Higgs states and the neu- 
tral pseudo Goldstone Higgs state depend on the mechanism for realizing 
small first family masses. If the r\ Higgs is the primary contributor to second 
family masses, the pseudo Goldstone and r\ Higgs states must have masses 
above 220 GeV, with numerical fits suggesting masses above 1 TeV, while 
if the r] Higgs is responsible solely for first family masses, the corresponding 
mass bounds drop to the range detectable at the LHC. We show that natu- 
ralness of small first family masses favors the latter alternative, and give an 
illustrative mass matrix texture model. 
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I. Introduction 

In a recent paper [1] we constructed extensions of the standard model, based on the 
hypothesis that the Higgs bosons also exhibit a threefold family structure, and that the flavor 
weak eigenstates are distinguished by a discrete Z 6 chiral symmetry that is spontaneously 
broken by the Higgs sector. Two models were analyzed in [1], the first with one three-family 
set of Higgs doublets, and the second with two three-family sets of Higgs doublets. In the 
three-Higgs doublet model, the leading cyclically symmetric approximation to the quark 
and lepton mass matrices has the "democratic" form with all matrix elements equal, leading 
to one massive and two massless fermion families. In the six-Higgs doublet model, for a 
wide range of Higgs potential parameters, CP is spontaneously broken, and this breaking 
simultaneously modifies the democratic Ansatz to give nonzero masses to an additional 
family (assumed in [1] to be the second family) in leading cyclic approximation. Corrections 
to the cyclic approximation were used in [1] to give first family masses, and a nontrivial 
CKM matrix. 

In performing numerical fits to the data using the models of [1], we took into account 
bounds on flavor changing neutral currents solely through the constraint provided by the 
Ki — K 2 mass difference, which led to strong asymmetries in the fits between the up and 
down quark sectors. M. Peskin [2] has pointed out the importance of including in the analysis 
experimental bounds on the D\ — D 2 mass difference, which is the up quark sector analog 
of the Ki — K 2 mass difference constraint. The purpose of this paper is to give the formulas 
and numerical results needed for this extension of the analysis of [1]. We also consider an 
alternative version of the model of [1], in which cyclic asymmetries in the <fi Higgs couplings 
are responsible for second family masses, and the r\ Higgs contributes significantly only to 
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first family masses. We show that this alternative is favored by requiring naturalness of 
small first family masses. In addition, we make the following three corrections to the model 
as originally formulated: (i) we correct the form of the CKM matrix, as pointed out in an 
Erratum [3] to [1], (ii) we include rephasings needed to make the diagonalized quark mass 
matrices positive real, and (iii) we correct combinatoric factors in the flavor changing neutral 
current amplitude (amounting to an overall factor of 2), and give a more accurate treatment 
of the hadronic matrix elements appearing in the flavor changing neutral current constraints. 

This paper is organized as follows. In Sec. II we give a synopsis of results needed 
from [1], including the corrections (i) and (ii) noted above. In Sees. Ill and IV we analyze 
K\ — K 2 and D\ — D 2 mixing induced by Higgs exchange, including the corrections (iii). 
In Sec. Ill we give formulas for calculating the Higgs, the rj Higgs, and the pseudo 
Goldstone Higgs contributions, in the six-Higgs doublet model, to both the K x — K 2 and 
D 1 — D 2 mass differences. In Sec. IV we use the formulas of Sees. II and III to derive 
a series of bounds on the Higgs masses, which are evaluated numerically using lattice and 
model calculations of the relevant hadronic matrix elements, for two possible mechanisms 
for realizing the first and second family masses. Irrespective of this choice of mechanism, we 
find that the <f> Higgs masses must be greater than 17 TeV, in accord with analyses [4] of 
generic multi-Higgs models. This bound also extends, by specialization to the case in which 
the rj Higgs couplings vanish, to the CP conserving case of the three-Higgs doublet model. In 
Sec. V we analyze the implications of requiring that small first family masses arise naturally, 
as opposed to arising by detailed cancellations between physically unrelated quantities, and 
give a simple mass matrix texture model corresponding to the case in which Yukawa coupling 
asymmetries are responsible for second family masses. In Sec. VI we repeat the numerical 
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fits of [1], taking into account the results derived in the preceding sections, and summarize 
our conclusions. 



II. Synopsis of Needed Results from the Six-Higgs Model 



The six-Higgs doublet model of [1] is based on the assumption that there are two 
discrete chiral families of Higgs bosons, 0„ and rj n , n — 1, 2, 3. These are coupled to discrete 
chiral families of fermions to give a Lagrangian that is exactly discrete chiral invariant, and 
that is approximately invariant under cyclic permutation of the discrete chiral components. 
The model is constructed so that the Higgs fields develop nonvanishing vacuum expectations 
in a CP violating phase, and it is assumed that the Higgs bosons couple much more strongly 
to fermions than the rj Higgs bosons, and similarly for their corresponding expectations 
(denoted respectively by and Q v ). As a zeroth order approximation to the model, only 
the Higgs expectations are retained and cyclic permutation symmetry is assumed, leading 
to a "democratic" mass matrix with one massive and two massless families, and a CKM 
matrix of unity. 

Deviations from cyclic symmetry, and the r) Higgs expectations, are then added back 
as a perturbation, giving as the Lagrangian mass term 




(la) 



f=u,d,e 



Here = diag(0, 0, 1) is the projector on the third family, is the Higgs Yukawa 
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coupling for flavor /, and a? is a 3 x 3 matrix of perturbations, given explicitly by 

°22 = \^22 + 3 Rf + + ^2 + ^<jf , 

(16) 

<T L=^L> w = exp(2?n/3) , cJ-cj* . 

In Eq. (lb) the terms /i£ m arise from small deviations from cyclic symmetry in the Yukawa 
couplings of the Higgs bosons, the terms 5^ arise from deviations from cyclic symmetry in 
the <p Higgs expectations, and the term W arises from contributions to the mass matrix of 
the weakly coupled 77 Higgs expectations. Further details of the structure of a* are given in 
Eqs. (38b) through (39c) of [1] and are used in Sees. V and VI below. However, the only 
property needed for the analytic calculations of Sec. Ill is that, since the Yukawa coupling 
asymmetries are all real because the Lagrangian in the six-Higgs doublet model is assumed 
to be CP invariant, the first order perturbations a* obey 

a{ 2 = a[{ . (lc) 

This restriction holds even though the model of [1] chooses a CP violating ground state. 
Defining 

M' f = 3M (3) + a f , (2a) 

we then [1] construct the bi-unitary transformation matrices U[ , for which U[M'j:U{i is 
diagonal, with the eigenvalues ordered in absolute value. The fermion mass eigenstate basis, 
up to rephasings to be discussed, is related to the primed basis by 

f L =u f L ] fr s , 

(26) 

f' R =u£fr s , / = M,e , 



and the CKM matrix Uckm is given [3] by 

C/ckm = • (2c) 

Since Eq. (2a) defines a degenerate perturbation problem, the matrices U[ R are 
constructed [1] in two stages: first the 2x2 submatrix of spanned by the first two families 
is diagonalized exactly, and then the solution to this problem is used to perturbatively 
construct the full 3x3 diagonalizing matrices. Because the analysis of flavor changing 
neutral current effects in the next section ignores third family mixings, it suffices for this 
analysis to discuss only the 2x2 submatrix diagonalization problem. Suppressing for the 
time being the flavor superscript /, we define the 2x2 submatrix m by 

m=( an au ) , (3a) 

which is brought to diagonal (but not necessarily real) form by matrices Vl,r, 

^=(o I) • (») 

with \k,i\ < 1^2 1- An explicit construction of V^r is given in Appendix B of [1]; the results 
obtained there can be simplified by using the symmetry of Eq. (lc) above, which (in terms 
of the quantities defined in Appendix B of [1]) implies that 





=A R = 


Wn 


| 2 + 


l^r I 2 


= A 


B L 


=B R = 


\&22 


| 2 + 


\sr I 2 

\ a 12\ 


= B 



(4a) 

ZL =(<7u + Cr 2 2)crl2 , Z R = (°"ll + cr 22) cr 12 > 
\zl\ =\Z R \ = \z\ . 

These relations, together with the results in Appendix B of [1], imply that 

_ / cos0 -exp(-# Lii? )sin0\ . . 

Vl ' r ~ Vexp(^ LiJ? )sin0 cos0 J ' [W) 
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with 

exp(i0 L>i? ) = z L:R /\z\ , (4c) 

and 

e = 5 ta »- 1 (j^) ■ (4rf » 

Although the construction just given suffices for the computation of the magnitudes 
of the CKM matrix elements, the calculation of the Higgs exchange amplitude in the next 
section requires care in the choices of phases. When we rephase the physical mass eigenstates 
/™i? s > the matrices V[ R transform according to 

Vl R ->Vl R = Dl R Vl R , (5a) 

with Z){ r diagonal unitary matrices. A correct choice of phases requires that the diagonal- 
ized mass matrix 

V[m f V^ (56) 

be real and positive; in other words, the matrices L>{ r must be chosen to absorb the phases 
of the diagonal matrix elements k[ 2 on the right of Eq. (3b) above. In addition, restricting 
ourselves now to the up and down quark flavor sectors, we shall require that the phase choices 
for the physical states put the matrix V^V^' , which is the 2x2 submatrix of the rephased 
CKM matrix when third family mixings are neglected, into the standard real form 

( Cl2 Sl2 ) , (5c) 

\-Sl2 C12 J 

with S12 and Ci 2 both nonnegative. These two phase requirements together fix the rephasing 
matrices D^ R up to an irrelevant overall phase. 
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To carry this construction out explicitly, we write Eq. (3b), in the up and 
quark sectors, as 

V f m f v n _ ( \4\ exp(^f) \ * — 7 

VLmVR ~\ |4|exp(^)J ' ] ~ UA ' 

and we write the adjoint of the 2x2 CKM matrix computed before rephasing as 



ydyu\ = ( ci 2 exp(i0n) s 12 exp(i9 12 )\ 
L L \ Sl2ex P(^ 21 ) Ci 2 exp(i0 22 ) ) 



with unitarity imposing the conditions 

C 12 + S 12 = 1 ) 

On + #22 = #12 + #21 (mod 2tt) . 
Then a simple calculation gives 

£>£ =diag[exp(-i0 n ),exp(-i0 21 )] , 

-Dj? =diag[exp« - i0 n ), exp(t6 d 2 - i0 21 )] , 

D U L =diag[l, exp(-i0 21 + i0 22 )] , 

D U R =diag[exp«), exp(i0£ - i0 21 + *0 22 )] . 
Corresponding to these, we find from Eqs. (4b) and (5a) that 



d= f exp(-i0n)cd -exp(-2<^ - i0 n )s d 
I exp(i(j)i - i6 2 i)s d exp(-i6 2 i)c d 



yd = f exp(-iO u + i6f)c d -exp(-i(f) d R -i9ii+i6f)s d \ 
Vr \exp{ict> R - %e 2X + 20 2 d )s d exp(-i0 21 + i6 d 2 )c d ) 

yu = ( c u -exp(-i(p u L )s u 

L \ exp(i(f>l - i9 2 i + i0 22 )s u exp(-?0 21 + i0 22 )c u 

V u = ( exp(i9^)c u -exp(-i(p R + i9^)s u 

R \ exp(i(f) R - i6 2 i + i0 22 + i9^)s u exp(-i0 2 i + i0 22 + i6^)c v 

with s d>u defined in terms of the angle 6 of Eq. (4d) by 

Cd,„ = cos(6 diM ) , s d , a = sin(e din ) . 



Equations (4a-d), (6a, b), and (7a-e) provide our starting point for calculating the Higgs 
exchange contributions to flavor changing neutral current processes. 



III. Analysis of Ki — K 2 and Di — D 2 Mixing 
Induced by Higgs Exchange 

We begin by extending the formulas of [1] for the Higgs exchange contribution to 
the K 1 — K 2 mass difference to the case when the rj Higgs and pseudo Goldstone Higgs 
contributions are also included, using the rephased matrices V^r of Eqs. (7a-e). Our starting 
point is Eq. (45c) of [1] for the AS = 1 terms in the Higgs Lagrangian density, calculated 
to zeroth order in the perturbation a? , which when extended to include the 77 and pseudo 
Goldstone Higgs couplings reads 

4 A c r = E E E {defU4L + B^)s + sefU4L + B$ 2llb )d\ 

£_=<P,rip=± F=R,I 

+de ( ^L + B% 2l5 )s + seg(Ag 21 + B% ll5 )d (8a) 

+d€p G {A PG i 2 + BpcWJ^S + S€p G (A PG 2l + B PG21 ~f 5 )d . 

The corresponding formula for the effective Hamiltonian density H^ s s=2 for the AS = 2 
process s + s — > d + d is 

Ki= 2 = - \ E E E + B$ 12l5 ) s -^d(A$ 12 + B$ 12l5 ) s 

t=<t>, VP =±F=R,I M iF 

~\d{Af R12 + B% 2lb )s-^ m d{Af Rl2 + B% 2l ,)s (8b) 

A M r i R 



I-,, . „ 1 

L PG 



d(A PG i 2 + Bp G12 -f 5 )s— — d(A PG i2 + BpGWjr^s . 

L, 1V1 r. 



The eight Higgs squared masses appearing in Eq. (8b) that carry superscripts (±) are equal 
in pairs, 

M^=M^f->, M 2 ^ = M§ + \ i = ^ . (8c) 
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Although they (as well as ) are given in terms of Lagrangian parameters by Eq. (46b) 
and Table II of [1], we will not use these expressions, but rather will treat the Higgs masses 
that are independent, after taking account of Eq. (8c), directly as parameters to be bounded. 

The subscripts 12 (or 21) in Eq. (8a) indicate the row 1 to column 2 (or row 2 to 
column 1) matrix element of the corresponding 2x2 matrix expressions for the A and B 
coefficients, which we now give. Because CP invariance of the Lagrangian for the six-Higgs 
doublet model implies that the Yukawa couplings appearing in Eq. (8a) are real, the A 
coefficients appearing in Eq. (8a) are related to the B coefficients as follows, 

A (±) _, B (±) 

ASlJ^Bpo , (8d) 



a 



a 



A PG = -i-^B 

n AV =(nj + n;)s 



(3) 



Defining 



M« =(J I) , m£ 2 =(° J) , p 3 =(j ."J 



the matrices -B^/ are given by 



5g =^(V?P3Vf - V^PsVf^ 
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(+) _ . 



(96) 



li 



(9c) 



the matrices B^ r are given by 

<' = -B# =^^[exp(<6l)V^M« 2 V? - exp(-i6l)V-Mg 2 Vf ] , 
5« = -Btf =^[eMM)V£M£ 2 V£ + exp(-^)^M( 1 x ) 2 ^ t ] , 

(3) 

and the matrices and 5p G are given by 

*g =X<[ exp ^^^ M ^^ - exp(-<6l)VgM 2 ( 5 2 vf ] , 
Bpg = ~ ^i^9 d v [eMiO)V L d M^ + expH^M&V?] . 

In the above formulas, 9 is the overall phase rotation angle between the and r\ Higgs 

expectations introduced in Eq. (21) of [1]. 

Taking the K to K matrix element of Eq. (8b), we get 

(K\H% S=2 \K) = -\s d A {K\(dsf\K) - \s d B {K\{d lb sf\K) , (10a) 



(9d) 



with 



I ^ M 2 (f) i + M 2(3) ^2 

\ 2^ 2^ 2^ M 2(p) I + 2(3) M 2 



PG 12 



(106) 



The corresponding formulas for the Higgs exchange contribution to the D to D transition 
amplitude are obtained by replacing if by D, d by u, and s by c in the above formulas, 
and replacing the explicit factors of i by — i in Eqs. (9b, c, d), with the latter substitution 
reflecting the fact that the up sector Yukawa couplings involve the charge conjugates of the 
Higgs fields. 

Substituting now the explicit forms given in Eqs. (7a-e) for the matrices V^'£, we get 
formulas for the sums Sj^, that determine the Higgs exchange contribution to the K to 
K transition amplitude, and for the corresponding sums S\ and S% that contribute to the 
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D to D transition amplitude . With an eye to how these formulas will be used in Sec. IV, 
we write them as 



S d A =s 2 d c 2 d eMi$ d )Pl , 



S d B = -4c d eM^d)P d B , 



S u A =s 2 u c 2 u exp(t$ u )PX , 



(11a) 



S u B = -sy u exp(t<S> u )P% , 



with the positive real quantities P A ™ B given by 



P d =3 



d\2 



(si) 



d\2 



+ -M) 



=3 



d\2 



(si) 



d\2 



+ -M) 



Pa =3 



(9 



u\2 



u\2 



+ tW 



Pr =3 



(5 



u\2 



U\2 



+ tW 



sin 2 y d , cos 2 y d 




sin 2 F u cos 2 r M 



M 



2(+) 



0i? 



M; 



2(-) 



0/? 



M 



1 1 

+ 



2(+) 



M. 



2(-) 



U\2 



r^ii r)R 

cos 2 F M sin 2 F u 



J_ J_\ 1 



u\2 



co^Xd j^sin 2 ^ 
M 2 ^ Way Ml G 



sin 2 X, , cos 2 X d 



My 



2(3) 



+ 



O 2 M 2 

\L AV lVlp G 



cos 2 X u t ttl sin 2 X u 



M, 



2(3) 



+ 



<t> 

n 2 AV mi 



PG 



sin 2 X u , til cos 2 X u 



M, 



2(3) 



+ 



O 2 M 2 

\L AV IVlpQ 



(lib) 



The mixing and phase angles appearing in Eqs. (11a, b) are given in terms of the various 
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phase angles denned above by 

Y d =\{dt + o d 2 ) , 

Xd^ + \tfR-<&)-\{oi + el) , 

<S> d =e*-0*-20 11 + 20 21 + , 

(12a) 

K=i(^ + ^) , 

$ u =0 « _ « _ 2 # 22 + 26 21 - (r R + 4>l) ■ 
As a check on our phase conventions, we note that when the model is CP conserving, 
which implies [1] the additional condition a 22 = cr(i, then Eq. (4d) reduces to 6 = 7r/4, and 
the following relations hold (modulo 7r), 

6 d,u = _ d,u = aTg(j d,u ? 

0f = -arg(7* M -arg(7* M , 

^ tt = -arg(7* tt + arg(7* M , (126) 

2011 =" 26 22 = <j) U L - K , 

26 12 = - 26 21 = -r L -<p d L ■ 
When substituted into Eq. (12a), these relations imply the vanishing (modulo n) of the phases 

$,i and $„. Consequently, with the phase conventions used in this paper, the imaginary parts 

of the Higgs exchange contributions to K — K and D — D mixing are a direct measure of 

the CP violating contributions to these amplitudes. 

IV. Higgs Mass Bounds 

We proceed now to derive bounds on the Higgs masses in the six-Higgs doublet 
model. Let AM* S and AM^ bs be respectively the measured value of the K\ — K 2 mass 
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difference and the experimental upper bound on the D\ — D 2 mass difference. Since it is 

reasonable to expect these to set upper limits on possible Higgs contributions to these mass 

differences, given respectively by [5] 

\AM^F K J =M^\(K\H^= 2 \K)\ , 

(13a) 

\AM^%J=M D 1 \(D\H^= 2 \D}\ , 
we get the basic inequalities 

AM? >M^\(K\n% s=2 \K)\ , 

(136) 

AM°£ S >M D l \(D\H^ =2 \D)\ . 
These inequalities will be used in this section, both independently and in combination, to 

derive a number of useful bounds on the Higgs masses. 

We begin by rewriting Eqs. (10a, b) and (11a, b) so as to exhibit the features that 

play a role in our various inequalities. Let us define px and as the negatives of the ratios 

of the scalar to pseudoscalar matrix elements appearing in Eq. (10a) and in its D meson 

analog, 

(K\(ds) 2 \K) 

Pk - 



(D\{ucf\D) 
(D\{u l5 cf\D) 

According to calculations of p K and p D by the vacuum insertion method [6, 7] and the MIT 
bag model [6], they are positive and small (roughly of order 0.1 in magnitude). The ratio px 
can also be extracted from lattice calculations that have been performed [7] for kaon matrix 
elements, giving the result px = 0.30 ± 0.05, again of positive sign. Although a similar 
lattice calculation is not yet available for the D system, we will assume that this follows 
the same pattern as observed in the K system, and that po (as suggested by the vacuum 
saturation and bag model calculations) is positive. Substituting Eqs. (10a), (11a), and (14a) 
into Eq. (13b), our two basic inequalities now take the form 
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M x AM° bs >^\(K\(d l5 s) 2 \K)\slc 2 d \P d B +p K Pi\ , 

M D AM° bs >^m(u l5 c) 2 \D)\slcl\P^+p D P2\ , 

with Pg U + pk,dPa U both sums of positive terms. Introducing the definitions 

2M K AM$ S 



(146) 



-Ed 



S ? 2 c? 2 |(K|(rf 75S )2|K)| 

2M D AM£ bs 
: S ? 2 c? 2 |(D|(lZ75c)2|D>| 

we rewrite the inequalities of Eq. (14b) as 



E K \ 2 > \sdCd\ 



(14c) 



1 

-Ed \ 2 > |s«c u | 



(Ud) 



\P%+PdP1\J "S12C12 
Although Eqs. (14d) are relevant for the numerical fits of Sec. VI, where the products 

\sdCd\ and \s u c u \ are known, they cannot be used to give fit-independent bounds on the Higgs 
masses, because either \s d c d \ or \s u c u \ can vanish. However, we shall now show that the 
sum \s d c d \ + \s u c u \ is bounded below by CKM matrix elements, permitting us to extract a 
useful inequality by combining the K meson and D meson flavor changing neutral current 
constraints. To see this, we substitute Eq. (4b) for V d,u , together with the definitions of 
Eq. (7e), into Eq. (6b) for the adjoint of the unrephased CKM matrix, and take absolute 
values of the matrix elements on the first row, giving 

S12 =\c u s d exp(-i(f) d L ) - s u c d exp(-i(f)l)\ < \c u s d \ + \s u c d \ , 

(15a) 

C12 =\c u c d + exp(i(p u L - i(p d L )s u s d \ < \c u c d \ + \s u s d \ . 
Multiplying these inequalities, we get 

S12C12 < \c d Sd\(s 2 u + c 2 u ) + \c u s u \(s 2 d + c 2 d ) = \c d s d \ + \c u s u \ , (156) 
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giving the needed lower bound. Adding the two inequalities in Eq. (14d), and using Eq. (15b), 
we get the master inequality 



Ek 



\P£+p K Pj\ 



E 



D 



\Pb+PdPa\ 



> 1 



(16) 



Since all terms in the denominators \P^ U + Pk,dPa U \ are positive, deleting any of 
these terms serves to make the left hand side of Eq. (16) larger, giving a number of simpler 
subsidiary inequalities that are consequences of the master inequality. Specifically, if we 
delete all terms in both denominators that do not refer to a given Higgs mass (i.e., if we set 
all of the other Higgs masses equal to infinity), we get a lower bound for the Higgs mass that 
we have retained; performing this in succession for the six Higgs masses we get the following 
inequalities, 



M (+) > 



M PG > 




E 



i) 



3(g d 4> ) 2 (sm 2 Y d + p K cos 2 Y d 



J K 



3(g d v ) 2 (l+p K ) 
2E 



+ 



4E 



D 



3(g") 2 (cos 2 Y u +p D sm 2 Y u ) 
E D 

3(gW (sin 2 Y u + p D cos* Y u ) 



n-l 

2 




3(^) 2 (1+Pd) 



K 



3«) 2 (sin 2 X, + ^cos 2 X,) 

2EkQ av 
3(g*ynl(cos 2 X d + p K sm 2 X d ) 



+ 



2E 



D 



I \ - 1 

2 



3(^) 2 (sin 2 X u +p D cos 2 X u ) 



2E D Q\ V 
3(g%) 2 {l 2 ( C os 2 X u + p D sm 2 X u ) 



2 



(17a) 



Applying the same procedure of successive deletion of denominator terms to the inequalities 
of Eq. (14d), we get a set of analogous inequalities [which, by use of Eq. (15b), imply those 
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of Eq. (17a)] that will be used in the numerical work of Sec. VI, 



M$ > max { r d 



> max < r d 



M$ > max { r d 



> max { r d 



M PG > max < r d 




3(^) 2 (sin 2 F d + p x cos 2 F d ) 
AE K 



2E 



AE 



D 



K 



3(^) 2 (sin 2 X d + p x cos 2 X d ) 



? '"it 



IE 



D 



3(^) 2 (sin 2 X u +p D cos 2 X u ) 



2E K tt 2 



AV 



3(g^Q 2 ( C os^X d + p K sm 2 X d ) 



2E D tt 2 



AV 



3(^) 2 n 2 (cos 2 X u +p D sin 2 X u ) 



r d = , r u = 

S12C12 S12C12 



(176) 



The bounds in Eq. (17a) still depend on the mixing angles X dyU and Y dyU defined 
in Eq. (12a); a set of (necessarily weaker) bounds that do not depend on these angles is 



obtained by using the inequalities, valid for p < 1, 

cos 2 Z + psin 2 Z =(1 — p) cos 2 Z + p > p , 

sin 2 Z + pcos 2 Z =(1 — p) sin 2 Z + p > p , 
giving (for Pk,d < 1) the inequalities 



< > 



(18a) 



M PG > 



E K 


1 
2 


E D 


M 1 

) 


3(9f) 2 PK 


+ 




) 2 Pd 




AE K 




1 
2 




AE D 


L3(^) 2 (1+^)J 


+ 


Wg u v 


) 2 (1+Pd) 



2E 



K 



Kg$) 2 PK 



+ 



2E 



D 



Hgfi 2 p D 



-1 



(186) 



2EkQ\v 

Kg^lPK 



+ 



2E D Vt 2 AV 

Hg^lPD 



1 \ —i 

2 
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Equations (17a, b) and (18b) are our final bounds for the Higgs masses in the six-Higgs 
doublet model. 

To obtain numerical values from the bounds of Eq. (18b), we first need to evaluate 
the ratios Ek,d defined in Eq. (14c). For Ek we use the measured value [8] AMj^ bs = 3.49 x 
1(T 12 MeV, together with M K = 497.7 MeV, M d = 6 MeV, M s = 115 MeV, f K = 160 MeV, 
C12 = 0.975, S12 = 0.221, and the lattice evaluation [7] 

\(K\(d l5S y\K)\ = 1.58 M 2 K f K , (19a) 

to give 

E K = 0.44 x 10~ 12 (GeV)- 2 . (196) 

For Ed, we use the experimental upper bound [8] AM D hs = 1.58 x 10 -10 MeV, together 
with M D = 1865 MeV, M u = 3.25 MeV, M c = 1.25 GeV, the lattice calculation [9] value 
Jd — 1-2/x, and the vacuum saturation approximation formula [7] 

\(K\(d l5 s) 2 \K) _ P K M 2 K[^Ml{M s + M d )-* + l] 

\(D\(u l5C y\D}\ / 2 M 2 [llM 2 (M c + M u )-2 + l] ' 1 ] 

to give 

E D = 27 x 10~ 12 (GeV)- 2 . (19d) 

Since the scalar to pseudoscalar ratio pd has not yet been computed on the lattice, we will 
assume that Pd = Pk = 0.3 in evaluating Eqs. (18b). 

To complete the computation of Higgs mass bounds, we need the values of the various 
Yukawa couplings appearing in Eqs. (18b). Here some assumptions about how the first and 
second family masses are generated are needed. If, following [1] , we assume that the Higgs 
expectations generate the third family masses, the 77 Higgs expectations generate the second 
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family masses, while cyclic asymmetries in the Yukawa couplings are responsible for the first 
family masses, then we get from Eq. (32b) of [1] the formulas 



M t d M b 



(2Ua) 



9 V on ' 9r) 



If we now assume equal Vt^ and Q v , so that = Q v = 71 GeV, we get the numerical values 

$=0.82, 4 = 0.020 , 

(206) 

$ =0.0059 , g$ = 0.00054 , 
giving the Higgs mass lower bounds 

M$>24TeV , 
M$ >470 GeV , 

(20c) 

>320 GeV , 
M PG >220 GeV . 

An alternative possibility, discussed in the next section, is that the second family 
masses are generated by cyclic asymmetries in the Yukawa couplings, with the first family 
masses generated by the 77 Higgs expectations. In this case, the second line of Eq. (20a) is 
replaced by 

<"l«V ^TW, • (21a) 

giving (for Q v = Q^) the numerical values ^ = 3.1 x 10~ 5 , g't = 5.7 x 10~ 5 , which imply 
the much weaker r\ and pseudo Goldstone Higgs mass bounds 

Mg } >5.5 GeV , 

M$ >3.7 GeV , (216) 

M PG >2.6 GeV . 
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Finally, we note that the first inequality of Eq. (18b) also applies to the CP conserving 
case of the three-Higgs doublet model of [1], for which = y/2 x 71 GeV, so that g^' u are 
a factor of y/2 smaller than given in Eq. (20b). This reduces the corresponding bounds of 
Eq. (18b) by a factor of y/2, giving for the CP-conserving three-Higgs doublet model the 
Higgs mass lower bounds 

Mg>17TeV . (22) 

In this model there are no rj Higgs states, and hence no possibility of neutral Higgs states 
that are not supermassive. 

V. Implications of Requiring Naturally Small First Family Masses 

We saw in the preceding section that the bounds on the rj and pseudo Goldstone 
Higgs masses depend on the coupling pattern assumed for the 77 Higgs discrete chiral triplet. 
In this section we classify possible r\ Higgs coupling patterns, based on a criterion of requiring 
naturally small first family masses. Referring to Eqs. (la) and (lb), we see that contributions 
to the mass matrix in the six-Higgs doublet model are of three distinct types, arising from 
deviations from cyclic symmetry in the <f> Higgs Yukawa couplings, deviations from cyclic 
symmetry in the <f> Higgs expectations, and contributions from the weakly coupled rj Higgs 
expectations. Since these three contributions are not directly related physically, detailed 
cancellations between them in the determination of the first family masses are a priori 
unlikely. Hence as a necessary (but not sufficient) condition for naturally small first family 
masses, we impose the condition that only one of these three contributions dominates in a 
leading approximation in which the first family masses are exactly zero. 

We begin by noting that the deviations 5[ from cyclic symmetry in the <fi Higgs 
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Yukawa couplings cannot dominate and lead to naturally zero first family masses. Let us 
suppose that the 8^ do dominate, and consider first the case in which the model chooses 
a CP conserving ground state, for which the parameters 8^ are all real. In this case the 
magnitudes 

\a{ x \ ~\8{ + u8{ + u)8{\ , 

(23) 

l°22l -|<*3 +W8{ + W8{\ , 

are equal, and so the first and second family masses are equal. Turning on a CP violation 
results in complex 8^s, for which the first and second family masses are no longer the same, 
but clearly a fine tuning of the amount of CP violation would be needed to achieve zero first 
family masses. Hence dominance of the Higgs expectation asymmetries 8^ is not compatible 
with naturally small first family masses. 

We consider next the case in which the contribution Rf of the rj Higgs expectations 
dominates, which is the scenario assumed in [1]. In this case the 77 Higgs expectations give 
rise to the second family masses, and the leading approximation to the first family masses is 
automatically zero, satisfying our criterion for naturally small first family masses. However, 
a potential problem arises when we examine the structure of the CKM matrix. In the leading 
approximation in which only Rf is retained in the mass matrix, the CKM matrix is unity. 
To get a nontrivial CKM matrix, we must add back the small perturbations 5 1 and [i\ m in 
Eq. (lb). According to Eqs. (4a-d), in each flavor channel we then get 

A-B~-\(T22\ 2 , 

|s| ~|<7 12 ||(7 22 | , (24a) 

which by the hypothesis of dominance of Rf is much less than unity. Hence, in particular, 
the up channel quantity s u = sin 6" is much less than unity in magnitude. But referring 
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now to the corrected expressions [3] for the CKM elements s± 3 and s 2 3, we have 

S13=|S3-<W3, S23 = |s 3 + 4|/3 , 

(246) 

S3 =C U « - < 3 ) , <*3 = S « eX P ( ) (^23 - ^23 ) , 

indicating that the spread of s 2 3 and S13 from their geometric mean is suppressed by the 
small quantity s u . This in turn requires relatively large parameter values <t 23 and/or a 13 
to give a satisfactory fit to the data, contradicting the starting assumption of a dominant 
Rf. We shall see evidence for this phenomenon in the next section, where we find Yukawa 
asymmetries comparable in magnitude to Rf , and hence substantial fine tuning in achieving 
small first family masses. 

We turn finally to the third case, in which the dominant contributions to the mass 
matrix come from the asymmetries /i{ m of the <p Higgs Yukawa couplings. For the leading 
approximation to the 2x2 submatrix m of the mass matrix, we then have (suppressing the 
flavor index /) 

3 VA*21 A*22 J ' V ^ 

with only two of the matrix elements in Eq. (25a) independent, since CP invariance of the 
Higgs Yukawa couplings implies [1] that 

Afci = M12 > ^22 = IHi ■ (256) 
In order for Eq. (25a) to have a zero eigenvalue, we must impose the additional condition 

I/U12I = |a*ix I H12 = exp(-ix)l4i > ( 26a ) 

an explanation for which must be sought in higher energy physics determining the Yukawa 
couplings. Taken together, Eqs. (26a) and (25b) imply that the matrix m takes the rank one 
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form 

3 \exp(ix)liii J 

with eigenvalues \ki\ =0 and | /«2 j = ||/in|, corresponding respectively to the first and second 
family mass eigenstates. ^From Eq. (4a), we find that the diagonalizing matrices V^r are 
given by Eq. (4b), with 

71 1 

9 = — , cos = sin 6 = — — , 

4 V2 (27a) 

<Pl =X , <Pr = X + 2 arg/in . 
Referring to Eq. (15a), we see that the sine of the Cabibbo angle si 2 is given now by 

Sl2 = i|expH0l)-exp(-^)| = |sini( X M - X d )| . (276) 

Averaging Eq. (24b) and the analogous expression obtained from the lower left corner of the 
CKM matrix, we get the following leading order expressions for si 3 and s 23 , 
Sl3 =|«3 - 41/3 , s 23 = \s 3 + d 3 \/3 , 

1 111 ( 27c ) 

S3=^K 3 -^3), d 3 = —eM-^X u + X d )}cos[-( X u -X d )M 3 -a u 23 ) , 

in which the coefficient of a 23 — <J 23 in d 3 is not now a small parameter. 

To complete the analysis of the third case, let us calculate the first family mass 
eigenvalue. There are four possible contributions to a nonzero first family mass: (i) deviations 
from the rank one condition of Eq. (26a) on the 2x2 submatrix m of the mass matrix, (ii) 
asymmetries in the (f> Higgs expectations S n , (iii) couplings to the third family through the 
mass matrix elements a 13 , a 31 and o" 2 3, o" 32 , and (iv) effects of the r] Higgs expectation 
term R in Eq. (lb). (We are continuing to suppress the flavor index / when not needed.) 
The simplest way to calculate the first family mass matrix eigenvalue \ki\ is to evaluate the 
absolute value of the determinant of the 3x3 mass matrix M' of Eq. (2a), which yields \k\\ 
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when divided by the product of the other two eigenvalues, giving to leading order in small 
quantities, 

M Jdet(3M(3» +ff )| 
3|/c 2 | 2 1 /in I 

To illustrate this in a simple texture model corresponding to the third case, let us assume 
that the contributions (i) and (ii) above are absent, that is, we assume that the condition 
of Eq. (26a) is exact and that cyclic asymmetries in the Higgs expectations are negligible. 
Then the only relevant contributions to the determinant in Eq. (28a) are (iii) and (iv) above, 
and a simple calculation gives the leading order formula 

3 1 

l«i| - 1 2^ ~ g[°"i3 _ ex P(#L)°"i3][o"3i - exp(-i<f> R )a 31 ]\ , (286) 

with 4>l,r given in Eq. (27a). When only the R term is retained, substituting Eq. (28b) into 
Eq. (la) yields the formulas of Eq. (21a). Within the simplified texture model, we have also 
calculated the CP violating angle o~i 3 appearing in the standard form [8] CKM matrix as a 
consequence of the CP violation carried by the R term. After considerable algebra, we find 

<f 13 ^ 2 v /2e u [cot(0^ - (j> d L ) - cot(0£ + 20i 3 )] - ^2e d csc(0£ - <j> d L ) , (29a) 

with the auxiliary quantities appearing in Eq. (29a) defined by 

n _ ^„i^d „u \ c u,d _ -9Re(/i 1 { R u4 ) 

U n - arg((7 13 - a 13 ) , e - — d — • (296) 

4 v 2 l^ii r 

To complete the specification of the texture model corresponding to the third case, 
we note that since the 2x2 diagonalizing submatrices V^r are maximally mixing in this 
case, wherever "1st" or "2nd" appears in the Higgs meson fermion family coupling Table II 
of [1], there now should appear "1st and 2nd", indicating couplings of equal magnitude of 
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first family to first family, first family to second family, and second family to second family. 
We also remark that the rank one condition of Eq. (26a) can be reexpressed as a model for 



the Yukawa asymmetries (3e m 


, by using the inversion formulas 




2 

= -Re(/in +/ii2 +/il 3 +/i3l) , 

y 


Pl2 


2 

= -Re(u;/iii + a7/ii2 + ^13 + ^31) , 
y 


Pl3 


2 

= -Re(a7/in + ufi 12 + /ii 3 + a7/i 3 i) , 


1321 


2 

=-Re(w^n +w^i2 +w^i3 + A*si) , 


(322 


2 

= -Re(cj/in + /ii2 + ^13 + w/i 3 i) , 


&3 


2 

=-Re(//n + 07/13 + 07/31) , 


(hi 


2 

= -Re(5fytii + u;/ii2 + 07/13 + //31) , 


(3 3 2 


2 

= -Re(/in + u;/i 12 + 07/13 + 07/31) , 


033 


2 

= -Re(cJ/!n +/ii2 + W/ii 3 + Ufl 31 ) . 



VI. Some Illustrative Numerical Fits 

We give in this section illustrative numerical fits, obtained by the method of min- 
imizing a "cost function" described in detail in Sec. IX of [1]. The mass and CKM cost 
functions are as in Eqs. (51b) and (52b) of [1], with the changes that we omit the term 
(f^ — il v ) 2 , which had little effect on the fits, and take the target values and standard devia- 
tions from the latest Particle Data Group [8] compilation. For the parameter cost function, 
we use Eq. (53a) of [1], with the changes that we omit the summation restrictions excluding 
the n = 3 and m = n = 3 terms, so that all asymmetries are treated symmetrically, and 
we take the exponent e as 2 rather than as 1. For the third case of the proceeding section, 
in which the second family masses arise from Yukawa coupling asymmetries, we also add to 
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the parameter cost function a term 




designed to keep the 77 Higgs contributions small, and start the iteration from preliminary 
parameter values determined with g? — 0, f — u, d, e. We omit flavor changing neutral cur- 
rent constraints from the cost function, so that there is no fine tuning to attempt to suppress 
flavor changing neutral current effects; instead we use the inequalities of Eq. (17b), evalu- 
ated using the parameters determined by the fitting procedure, to give lower bounds on the 
Higgs masses that guarantee sufficiently small flavor changing neutral current contributions 
to K — K and D — D mixing. 

Fitting results for the second case of the proceeding section are given in Table I, 
obtained with a standard deviation for the Yukawa asymmetries (5 of o"p arame ter = 0.02. For 
this fit, the maximum \(3\ values in the up, down, and electron sectors are 0.032, 0.058, 
and 0.010 respectively. Although the iteration is started with the 77 Higgs coupling values of 
Eq. (20b), the converged fit has significantly smaller couplings = 0.0031 and = 0.00039, 
indicating that the Yukawa asymmetries make substantial contributions to the second family 
masses. This means that the 77 Higgs contributions do not dominate the Yukawa asymmetries, 
and thus substantial fine tuning is involved in achieving small first family masses. For the 
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fit of Table I, the Higgs mass bounds obtained from Eq. (17b) are 

M$ >520 TeV , 

M ( ~J >440 TeV , 

M$ > 1.3 TeV , (32a) 
M$ >1.1 TeV , 
M PG >1.0 TeV . 

In Table II we give fitting results for the third case of the preceding section, obtained 
now with a standard deviation for the Yukawa asymmetries (5 of o" paramet er = 0.08. For this 
fit, the maximum \(3\ values in the up, down, and electron sectors are 0.042, 0.052, and 0.11 
respectively, with the relatively large Yukawa coupling asymmetry needed in the electron 
sector reflecting the fact that in this sector the ratio of the second family to third family 
mass is biggest. For the fit of Table II, the Higgs mass bounds obtained from Eq. (17b) are 

M$ >370 TeV , 

M$ >620 TeV , 

Mg } >210 GeV , (326) 
M$ >220 GeV , 
M PG >140 GeV . 

To conclude, in order for flavor changing neutral current effects in our models to 
be sufficiently small, the Higgs masses must be very large, far outside the regime in 
which conventional perturbative Higgs physics applies (see [10] for a recent review of both 
perturbative and strongly coupled Higgs models). Our results are consistent with general 
group theoretic analyses of flavor changing neutral currents in multi-Higgs doublet extensions 
of the standard model [11], which when applied to our models imply that flavor changing 
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neutral currents cannot cancel kinematically, but must be eliminated either by fine tuning 
(an option we have ruled out by the inequalities of Sec. IV) or by having some very large 
Higgs masses. From an experimental viewpoint, the most interesting scenario within our 
models corresponds to the third case discussed in Sec. V, in which the r\ Higgs couplings 
are small enough that their expectations play an important role only in determining the 
first family masses and in giving rise to CP violation. In this case the i] Higgs and pseudo 
Goldstone Higgs states are permitted by our mass bounds to be light enough to be seen in 
experiments at the LHC. The simplest model of this type would be one in which the and 
i] Higgs self-interactions have similar structures, with a weak — rj coupling. Massiveness of 
the Higgs states would then imply massiveness of the corresponding r\ Higgs states, with 
only one neutral and two charged pseudo Goldstone Higgs states potentially observable at 
LHC energies. In such models, one simultaneously has observable light Higgs states (the 
pseudo Goldstone triplet) and "new physics" implied by the strongly self-coupled and rj 
Higgs sectors. 
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Table I. Six-Higgs doublet model fit to experimental data: 
second case of Sec. V 



quantity target value fitted value 



M u 


0.0033 


0.0033 


M c 


1.25 


1.26 


M t 


173.8 


174.0 


M d 


0.0060 


0.0064 


M s 


0.115 


0.111 


M b 


4.25 


4.24 


M e 


0.00051 


0.00051 




0.1057 


0.1057 


M T 


1.777 


1.777 



S\2 

sin 5 13 \ 



0.221 
0.0059 
0.039 
0.60 



0.221 
0.0062 
0.037 
0.48 
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Table II. Six-Higgs doublet model fit to experimental data: 

third case of Sec. V 



quantity target value fitted value 



M u 


0.0033 


0.0032 


M c 


1.25 


1.25 


M t 


173.8 


173.9 


M d 


0.0060 


0.0065 


M s 


0.115 


0.096 


M b 


4.25 


4.25 


M e 


0.00051 


0.00051 




0.1057 


0.1057 


M T 


1.777 


1.777 


Sl2 


0.221 


0.221 


Sl3 


0.0059 


0.0059 




0.039 


0.039 


sin S 13 \ 


0.60 


0.33 
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